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Abstract 

Perturbed angular correlations were measured for Gd ions implanted into gadolin- 
ium foils following Coulomb excitation with 40 MeV 16 O beams. A technique for 
measuring the relative magnetizations of ferromagnetic gadolinium hosts under in- 
beam conditions is described and discussed. The combined electric-quadrupole and 
magnetic-dipole interaction is evaluated. The effect of nuclei implanted onto dam- 
aged or non-substitutional sites is assessed, as is the effect of misalignment between 
the internal hyperfine field and the external polarizing field. Thermal effects due to 
beam heating are discussed. 
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1 Introduction 



Gadolinium foils are used extensively for in-beam measurements of hyperfine 
interactions and nuclear moments [1,2,3,4,5,6,7,8,9,10,11,12,13]. Magnetized 
gadolinium foils are used in preference to iron foils in many applications of 
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the transient-field technique [5,13] to measure nuclear g factors because larger 
perturbations of the particle-7 correlations can be obtained under otherwise 
similar experimental conditions [6,7]. One disadvantage of gadolinium hosts, 
however, is that the magnetization is not well controlled, being sensitive to 
the crystalline structure of the foil, and varying considerably with both the 
applied field and the temperature, even at temperatures well below the Curie 
temperature of 293 K [14]. 

As cooling with liquid nitrogen (77 K) is convenient in accelerator laboratories, 
most in-beam measurements which employ gadolinium hosts have been per- 
formed at somewhat higher temperatures near 90 K because of thermal losses 
and the effects of beam heating on the target. In transient-field measurements 
an external magnetic field, with strength typically in the range from 0.05 to 
0.1 T, is applied to polarize the gadolinium foil. Higher fields are avoided to 
ensure that bending of the primary beam is negligible [15,13]. The spacial 
profile of the polarizing field along the beam direction is designed to minimize 
beam bending effects rather than to produce a uniform field at the target lo- 
cation. Depending upon the design of the pole pieces and the location of the 
target, the profile of the external field across the target may not be uniform. 
Under these conditions the magnetization in the beam spot must be carefully 
related to the off-line magnetization measurements. 

In this paper a method of determining the relative magnetization of gadolinium 
foils under in-beam conditions is described. The static hyperfine magnetic field, 
which acts on the nuclei of Coulomb excited 2f states in 1 54,i56,i58,i60Q c j ! j g uge( j 
to probe the local magnetization at the beam spot under in-beam conditions. 

Measurements similar to the present work were performed by Skaali et al. [1], 
and Kalish et al. [3] . Additional complementary information was also obtained 
by Hausser et al. [6]. In these previous works, however, the focus was on the 
precessions of the 4+ and 6f states as probes of the strength of the transient 
hyperfine magnetic field. Here the emphasis is on the 2f states and the strength 
of the static hyperfine magnetic field. 

While the interpretation of the hyperfine interactions in terms of the relative 
magnetization of different samples turns out to be rather straight-forward, 
there are several additional phenomena associated with the hyperfine fields 
and the ion-implantation process that may have bearing on the interpretation 
of the data. The present work will therefore include discussions of: 

(1) the presence of the electric- field gradient in the gadolinium matrix, which 
means that the hyperfine interaction cannot universally be treated as a 
pure magnetic interaction 

(2) the effects of those implanted nuclei which reside on damaged or other 
non-substitutional sites 
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(3) the magnitude of the transient hyperfine magnetic field, which acts on 
the implanted ions as they slow within the host, and which for Gd in 
gadolinium has the opposite sign to the static hyperfine field 

(4) a possible misalignment between the hyperfine magnetic field and the 
external polarizing field, which can affect the interpretation of the per- 
turbed angular correlation when the host is not fully saturated 

The paper is arranged as follows: The next section (section 2) reviews pre- 
vious work on the electric and magnetic hyperfine fields experienced by Gd 
ions in gadolinium. Section 3 describes the in-beam measurements, presents 
examples of 7-ray and particle spectra, and summarizes the excited-state life- 
times, g factors and quadrupole moments adopted for the analysis. The off-line 
magnetization measurements are presented in section 4. Section 5 concerns 
the 'unperturbed' angular correlations measured at room temperature. It in- 
cludes a review of the formalism, along with the results and a discussion of 
the measurements, which indicate the direction of the electric-field gradient in 
the gadolinium foils. The perturbed angular correlation results are presented 
in section 6. The formalism and data analysis procedures are described, the 
effects of combined electric-quadrupole and magnetic-dipole interactions are 
evaluated, as are the effects of nuclei on damaged sites. The results are dis- 
cussed in Section 7. 



2 Hyperfine fields in gadolinium hosts 

Gadolinium has a hexagonal close packed (hep) crystal structure. Nuclei within 
the gadolinium matrix therefore experience a hyperfine electric-field gradi- 
ent along with the magnetic dipole interaction. Below the Curie temperature 
of 293 K, gadolinium is a simple ferromagnet with a magnetic anisotropy 
that has a complex dependence on temperature. Thus the easy direction of 
magnetization chane; function of temperature, as does the magnetic 

hyperfine field. Mossbauer studies [2] reveal anisotropic magnetic hyperfine 
interactions in gadolinium crystals, with a field of \B st \ = 37.3(5) T at 4.2 K 
oriented at 28° to the c axis (unpolarized sample). The electric field gradi- 
ent is apparently less sensitive to temperature [8]. The measured splitting at 
4.2 K for 155 Gd, eqQ/h = 108 ± 1 MHz, implies an electric field gradient 
of V zz = 3.44(6) x 10 17 V/cm 2 , assuming Q = 1.30(2) b. This electric field 
gradient is consistent with that found for high-spin isomers in 147 ' 148 Gd at 
temperatures near 400 K [8]. 

The electric-quadrupole and magnetic-dipole hyperfine fields acting on Gd iso- 
topes in various gadolinium samples at different temperatures have been stud- 
ied previously by many techniques (see Refs. [2,4,8] and references therein). In 
the present work the Mossbauer values at 4.2 K will be taken as the point of 
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reference. The sign of the static magnetic field has been determined to be neg- 
ative, for example from in-beam perturbed angular correlation measurements 
[1,3,6]. 

Since the present work concerns in-beam implantation of Gd into gadolinium, 
the transient hyperfme magnetic field also acts on the ions before they come 
to rest. Thus the net perturbation of the nuclear spin distribution has con- 
tributions from both the transient and static hyperfme magnetic fields. For 
states of spin 4 + and higher in the ground-state band, the electric quadrupole 
interaction can be ignored (see below and [6]) and the perturbation of the 
angular correlation is manifested essentially as a rotation through the angle 

AQ = uT + A9 tU (1) 

where A9 t { is the precession angle due to the transient field and lot is that 
due to the static field. Furthermore, 

lot = -g — B st t, (2) 

and 

A6 ti = -g^- f Ts Btf(t) e-'/ T (ft, (3) 
n Jo 

where r is the meanlife, g the g factor of the excited nuclear state, and T s is 
the time taken for the recoiling ions to stop in the ferro magnet; B st and B tT are 
the static- and transient-field strengths, respectively. Since -B st and B tT have 
opposite signs for Gd in gadolinium, the two contributions tend to cancel. 

Equation (1) is correct only for small precession angles because the static-field 
perturbation includes an attenuation as well as the rotation of the radiation 
pattern. Furthermore, for the longer-lived 2f states the quadrupole interaction 
cannot safely be ignored. The formalism needed for a rigorous analysis of the 
data is presented in section 6 along with a discussion of the combined effect 
of the electric quadrupole and magnetic dipole interactions on the observed 
angular correlations. 



3 IMPAC Measurements 

3. 1 Experimental Procedures 

Hyperfme fields acting on Gd ions implanted into gadolinium were measured 
using the Implantation Perturbed Angular Correlation (IMPAC) technique, 
following procedures similar to those in an earlier study of Pt in gadolinium 
[11]. The measurements were performed using 40 MeV 16 4+ beams from 
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Table 1 



Summary of experiments. T is the nominal temperature on the target frame, /beam 
is the beam current. Target A is 16.9 mg/cm 2 thick. Target B is 6.2 mg/cm 2 thick. 



Run 


Target 


T 

(K) 


beam 

(unA) 




Type of measurement 


I 


A 


300 


2.5 


W(ftv): 9^ -- 


= 0°, ±15°, ±30.5°, ±45°, ±55°, ±65° 


II 


A 


90 


2.5 


Precession: 


e(±65°), e(±120°) 










W(9^): # 7 = 


= 0°, ±15°, ±30.5°, ±45°, ±55°, ±65° 


III 


A 


90 


0.375 


Precession: 


e(±65°), e(±120°) 


IV 


B 


90 


2.5 


Precession: 


e(±65°), e(±120°) 










W(9^): # 7 = 


= ±31°, ±45°, ±65° 



the Australian National University 14UD Pelletron accelerator. Table 1 gives 
a summary of the angular correlation and nuclear precession measurements 
performed. As will be discussed below, the temperature shown in Table 1 is the 
nominal temperature of the target frame, which does not necessarily represent 
the temperature at the beam spot. 



The two targets employed were the same as those used in a recent study of 
transient-field strengths for high-velocity Ne and Mg ions traversing gadolin- 
ium hosts [16]. Target A consisted of a rolled and annealed 16.9 mg/cm 2 
thick gadolinium foil. To aid with thermal conduction, the gadolinium foil 
was sandwiched between two 12 /xm thick indium-coated copper foils hav- 
ing 6 mm diameter holes punched through at the beam position. Target B 
consisted of 0.1 mg/cm 2 nat C, a thin flashing of copper (0.02 mg/cm 2 ) to as- 
sist adhesion, 6.2 mg/cm 2 of rolled and annealed gadolinium, and a 'thick' 
(5.65 mg/cm 2 ) copper backing. Both gadolinium foils were cold rolled, be- 
ginning with 0.025 mm thick foil of 99.9% purity purchased from Goodfellow 
Cambridge Limited. After rolling they were annealed in vacuum at ~ 800° C 
for ~ 20 min. To provide additional support and thermal conduction, target 
B was attached to a 12 jum thick copper foil using ~ 1 mg/cm 2 of indium 
as adhesive. The beam entered the carbon side of target B. For both targets 
the Coulomb-excited nuclei of the gadolinium layer recoiled with energies of 
~ 11 MeV and subsequently stopped within the gadolinium layer. 

Backscattered 16 ions were detected in two silicon photodiodes, masked to 
expose a rectangular area 8.5 mm wide by 10.2 mm high, and placed 23.7 mm 
from the target 4.0 mm above and below the beam axis at back angles. The 
backscattered particle spectrum extended from ~ 27 MeV, due to scattering 
at the front surface of the Gd foil, down to ~ MeV due to scattering in 
the depth of the target. The threshold was set at ~ 3 MeV. The 12 C layer of 
target B produced some a-particle groups from 16 O + 12 C reactions, which 
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allowed an in-beam calibration of the particle detector, but did not otherwise 
interfere with the measurement. The beam species and energy were chosen 
to ensure that 7 rays detected in coincidence with backscattered beam ions 
originate from Coulomb-excited Gd nuclei that recoil and stop well within the 
gadolinium layer. 

Gamma rays were detected using two ~ 20% efficient detectors placed ~ 7 cm 
from the target and two ~ 50% efficient high-purity Ge detectors placed 
15.2 cm from the target. The larger Ge detectors were placed at # 7 = ±120° 
to the beam axis throughout the measurements. For the precession measure- 
ments, the forward Ge detectors were placed at 9 1 = ±65°, near the maximum 
slope of the particle-7 angular correlation. These detectors were also moved 
through a sequence of angles to measure the angular correlations; see Table 1. 

3.2 Particle and 7 -ray spectra 

Particle and 7-ray spectra for the two targets are shown in Fig 1. The a 
particle group (s) at low energies, appears in the spectrum for target B due to 
reactions on the Carbon layer. The corresponding 7-ray spectra, measured in 
coincidence with the detected particles, are essentially identical in the region 
of interest, i.e. below 250 keV. 

Table 2 lists the observed 7-ray transitions having energies between 70 keV and 
250 keV. Relative 7-ray intensities at 6 1 = 65° are also given. The uncertainty 
in the relative intensities is better than 5% of the quoted value for the strongest 
lines and ~ 20% for the weakest. The 2f — > 0^ transitions in the even isotopes 
are sufficiently resolved from the close-by 7/2^ — > 5/2]" transitions in 155 Gd 



and 157 Gd. 



Table 2: 7-ray lines from Coulomb excitation of nat Gd 



with 40 MeV 16 0. 



£ 7 (keV) a Nucleus If -> Ij 



Intensity 6 



75.3 160 Gd 2+ -> 0+ 

76.9 157 Gd 7/2- -> 5/2 

79.5 158 Gd 2+ -> 0+ 

86.6 155 Gd 7/2- -> 5/2 
89.0 156 Gd 2+ -> 0+ 

95.7 157 Gd 9/2- -> 7/2 
105.3 155 Gd 9/2- -> 7/2 



57.9 



25.0 



80.4 



17.6 



100 



2.8 



2.6 



6 



Table 2: (Continued) 



E 1 (keV) a Nucleus If -> // Intensity h 

11/2" -> 9/2" 1.8 

2+ -> 0+ 22.0 

7/2~ -> 3/2- 5.3 

11/2" -> 9/2" 0.9 

7/2- -> 3/2- 7.4 

9/2" -> 5/2" 2.6 c 

4+ -> 2+ 33.0 

4+ -> 2+ 31.4 

9/2- -> 5/2- 2.5 

4+ -> 2+ 22.8 

11/2" -> 7/2" 1.4 

11/2- -> 7/2" 1.6 

4+ -> 2+ 1.6 

a 7-ray energies from Refs. [17,18,19,20,21,22]. 

b Relative 7-ray intensity at # 7 = 65°. 

c This transition was not directly observed; see text. 



120.1 


157 Gd 


123.1 


154 Gd 


131.4 


157 Gd 


140.6 


155 Gd 


146.1 


155 Gd 


172.8 


157 Gd 


173.2 


160 Gd 


181.9 


158 Gd 


191.7 


155 Gd 


199.2 


156 Gd 


215.6 


157 Gd 


246.2 


155 Gd 


247.9 


154 Gd 



The 172.8 keV, 9/2" — > 5/2" transition in 157 Gd was not observed, but its 
presence was inferred from the observation of the 9/2" — > 7/2~ transition and 
the known branching ratio [20]. A small correction to the precession data was 
made to account for the effect of the overlap between the 173 keV 4| — > 2+ 
transition in 160 Gd and this much weaker transition in 157 Gd. 



3. 3 Adopted lifetimes, quadrupole moments and g factors 



The level lifetimes, quadrupole moments and g factors required for the follow- 
ing analysis are summarized in Tables 3 and 4. The present analysis assumes 
that, in each isotope, g(2f) = g{^t) = 9(^1) ■ This assumption is justified 
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100 200 300 

energy [keV] 



Fig. 1. 7-ray spectra recorded in coincidence with backscattered beam ions (see 
insets) for the two gadolinium targets. Target B has a thin layer of Carbon which 
gives rise to the a-particle group. See Table 2 for detailed identification of lines 
between 70 and 250 keV. 

Table 3 



Lifetimes in the Gd isotopes 



Nucleus 


Qo a (b) 


level 




mean life (ps) 










rotor 


experiment 


Ref. 


154 Gd 


6.25 


2+ 


1710 


1710(20) 


[23] 


156Q d 


6.83 


2+ 


3270 


3270(60) 


[23] 


156Q d 


6.83 


4+ 


160 


161.5(2.5) 


[19] 


158 Gd 


7.10 


2+ 


3730 


3730(70) 


[23] 


158Q d 


7.10 


4+ 


219 


214(3) 


[21] 


160Q d 


7.27 


2+ 


3910 


3910(80) 


[23] 


160 Gd 


7.27 


4+ 


258 







a Intrinsic quadrupole moments Qo are from Ref. [23]. 



both by experimental evidence and theoretical expectations [9]. 

The necessary lifetimes and quadrupole moments have been measured for all 
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Table 4 



Quadrupole moments and g factors of 2f states in the Gd isotopes 



Nucleus 


E x (keV) 


i, 

rotor 


Q (b) 

experiment c 




154 Gd 


75 


-1.79 


-1.82 (4) 


0.430 (30) 


156 Gd 


79 


-1.95 


-1.93 (4) 


0.387 (4) 


158 Gd 


89 


-2.03 


-2.01 (4) 


0.381 (4) 


160 Gd 


123 


-2.08 


-2.08 (4) 


0.364 (17) 



a Adopted g factors from the tabulation in Ref . [24] . 

6 Rotor model spectroscopic quadrupole moments derived from the Qo values in Table 3. 
c Quadrupole moments are from muonic atom x-ray measurements [25]. 



but the 4+ state of Gd. Since the rotor model gives an excellent description 
of the measured lifetimes and quadrupole moments in 156 Gd and 158 Gd, and 
160 Gd is even more deformed than these isotopes, the rotor lifetime is adopted 
for the 4+ state in 160 Gd. Experimental values are used in all other cases. 



4 Magnetization measurements 



The magnetizations of the gadolinium foils were measured off line with the 
Rutgers magnetometer [26]. For these measurements samples were cut from the 
same rolled and annealed foils as those used to make the targets. The results 
are summarized in Table 5. For these foils, and many similar foils prepared by 
rolling and vacuum annealing, the magnetization is found to vary with both 
temperature and applied field. Within a few percent the magnetizations of the 
present foils track those of a single crystal magnetized along the b axis [14] . 

Rolled and annealed thin Gd foils, as used in nuclear experiments, typically 
have a texture such that the foil resembles a quasi-single crystal with the c 
axis perpendicular to the plane of the foil (i.e. the basal planes of the micro- 
crystals are in the plane of the foil) [27]. This texture is confirmed by x-ray 
diffraction measurements [28], and by the fact that the magnetization curve 
versus temperature resembles that of a single crystal magnetized along the 
b axis. Fig. 2 shows the results of a detailed study of the magnetization of a 
5.3 mg/cm 2 Gd foil as a function of temperature, which is very similar to that 
of a single crystal magnetized along the b axis [14] . 

The variation of the single-crystal magnetization with the external field is 
shown in Fig. 5 of Ref. [14] . In the region around 90 K a factor of two change 
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Table 5 



Results of off-line magnetization measurements a . 



O 1 

b ample 


lemp 


^ext 


M 


a 


0"sc 




(K) 


(Oe) 


(Gauss/cm 3 ) 


(Gauss/g) 


(Gauss/g) 


A 


100 


500 


1523 


193 


~ 205 


B 


50 


783 


1400 


177 


~ 190 


B 


77 


633 


1618 


205 


~ 200 



a M is the magnetic moment per cm 3 while a is the mag- 
netic moment per gram, thus M = ap, where p is the density 
in g/cm 3 . a sc is the magnetization of a single crystal mag- 
netized along the b axis as read off Fig. 5 in Ref. [14]. The 
quantities are given in cgs units for ease of comparison with 
Ref. [14]. The uncertainties in the measured magnetizations 
are of the order of 3%. There is an uncertainty of the order 
of 5% associated with reading a sc from the figure. 



0.20 




q Qr^ I i i i i I i i i i I i i i i I i i i i I i i i i 

50 100 150 200 250 
T (K) 



Fig. 2. Measured magnetization as a function of temperature for a 5.3 mg/cm 2 
gadolinium foil [12]. The line is drawn to guide the eye and the error bars represent 
the uncertainty in the absolute magnetization. The external polarizing field was 
0.09 T. This magnetization curve is very similar to that of a single crystal mag- 
netized along the b axis [14]. The arrow indicates the approximate temperature at 
which many in-beam measurements are performed. 



in the external field (from 0.05 to 0.1 T) changes the magnetization by only 
5%. The primary focus of this paper is therefore on the temperature variation 
of the magnetization. 
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Fig. 3. Left: Schematic of particle detector defining the co-ordinate frame. The beam 
is along the z-axis and the magnetic field direction is along the y-axis. The 7-ray 
detectors are in the xz-plane. Right: Definitions of the spherical polar angles (/3, 7) 
which specify the direction of the electric field gradient (EFG), V zz , with respect to 
the magnetic field direction and the beam axis. The EFG in the gadolinium target 
foils studied here is aligned predominantly along the beam axis, i.e. (5 = 90°, 7 = 0°. 

5 Unperturbed angular correlations 

5. 1 Formalism 

The theoretical expression for the unperturbed angular correlation after Coulomb 
excitation can be written as (see Ref. [29,30] and references therein) 

W(9 7 , 7 ) = £ y/(2k + l)( Pkq (8, 4>))F k Q k D k g *(^, 7 , 0), (4) 

k,q 

where {p kq {9,(j))) is the statistical tensor, which defines the spin alignment of 
the initial state, and which depends on the particle scattering angles (9, 0) and 
the geometry of the particle detector. F k represents the usual F-coefficients 
for the 7-ray transition, Q k is the attenuation factor for the finite size of the 
7-ray detector, and -Dg O (0 7 , 9 7 , 0) is the rotation matrix, which depends on the 
7-ray detection angles (0 7 ,</> 7 ). In the applications of interest k = 0,2,4. The 
co-ordinate frame is right-handed with the beam along the positive z-axis as 
shown in Fig. 3. In the present work the 7-ray detectors are in the xz plane, 
thus 7 = 0, and the rotation matrix is equivalent to an associated Legen- 
dre polynomial. The statistical tensors, and hence the unperturbed angular 
correlations, can be calculated accurately for the reaction geometry using the 
theory of Coulomb excitation. 

The Coulomb excitation calculations performed here are based on the de Boer- 
Winther code [31]. In this code the statistical tensors are evaluated in the 
particle-scattering plane. To calculate the angular correlations requires the 
tensors corresponding to scattering at angle 4> as defined in Fig. 3. These are 
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given by 



MM) = E/MM)£j,(&o,o) 



(5) 



= p kq (9,0)e' 



(6) 



where pk q (9,0) are from the de Boer-Winther calculation. Thus the required 
average statistical tensor at a given beam energy is given by 



where the integrals are over the dimensions of the particle detector and ^ is 
the cross section for Coulomb excitation corresponding to the scattering angle 
9. In the geometry used here (Fig. 3) there are two particle detectors placed 
symmetrically about the beam axis such that the numerical integration can 
be limited to the positive quadrant, 0° < <fi < 90°. The factor e iq ^ can then be 
replaced by (e <9 * + e~ iq4> + e^^ +7r) + e - i<i,w>+7r) )/4, which is cosg0 if q is even 
and is zero if q is odd. 

To obtain the statistical tensors of direct relevance to the present experiments, 
a further integration was performed to average over the energy loss of the beam 
in the target. A correction for feeding from higher states in the ground-state 
band was also made. Since the feeding path is only along the ground-state 
band, the statistical tensor of the fed state i can be evaluated iteratively using 



where pk q {i) is the unfed statistical tensor for the state and p dircct (^) j s the 
direct population of the state by Coulomb excitation. p tota1 ^ — |— 1 ) is the total 
population of the ground-band level above level i, including direct excitation 
and feeding contributions, if any. U k {i + 1 — > i) is the [/-coefficient for the 
% + 1 — > % transition [32]. These feeding corrections are small in the present 
work. In all cases feeding from states with / > 2 in the ground-state band 
contributes less than 7% of the total intensity in the 1\ — > 0^ transition. 

The resultant nonzero statistical tensors for the 2f state of 156 Gd are shown 
in Table 6. For an annular counter only the q = tensors are non zero. The 
broken azimuthal symmetry gives rise to the finite pt q values for q ^ 0, however 
these terms are small in the present case because the scattering angle remains 
near 180° and the spin of the excited state is aligned predominantly in the 
plane perpendicular to the beam. 



{Pkq) 



(7) 




Pk q (i)P dircct (i) + Uk(i + 1 -* i)pk q (i + l)P total (i + 1) 

^direct _|_ ^3total^ _|_ 



(8) 
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Table 6 

Statistical tensors for the 2f state of 156 Gd. 



k 


9 


Pk,q 








1.0000 


2 


2 


0.0061 


2 





-0.5033 


4 


4 


0.0002 


4 


2 


-0.0304 


4 





0.4404 



5.2 Angular correlation results and discussion 



Figures 4 and 5 show examples of comparisons between the calculated un- 
perturbed angular correlations and the data. The angular correlations for the 
2\ — > 0^ transitions shown in Fig. 4 were obtained at room temperature (run 
I), above the Curie temperature. Thus there is no magnetic-dipole perturba- 
tion, but there is expected to be an electric field gradient (EFG). In Fig. 4 
the dotted lines indicate the angular correlation anticipated for a target with 
the expected electric field gradient of V zz = 3.44 x 10 17 V/cm 2 distributed 
isotropically. There is no evidence for any attenuation of the anisotropy due 
to electric field gradients. This observation again confirms that these rolled 
and annealed foils have a texture such that the basal planes are perpendicular 
to the beam, in the plane of the foil. There is no perturbation because the c 
axis, and hence the EFG, is along the beam direction, perpendicular to the 
plane of spin alignment. 



The data shown for the 4+ — > 2± transitions in Fig. 5 include data obtained 
at both room temperature and at 90 K (runs I and II), where magnetic per- 
turbations are present. The calculated and measured angular correlations are 
in agreement. There is no observable difference between the measurements at 
90 K and 300 K because the unperturbed angular correlations can be recov- 
ered by adding together the data for both directions of the external field. This 
procedure cancels out the relatively small rotation of the radiation pattern 
which changes direction when the external field direction is reversed. (The 
same procedure cannot be applied for the 2f states because the perturbations 
are much too large.) 
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2.0 



1.5 - 



1.0 



0.5 



0.0 
2.0 



1.5 - 



1.0 



- 1 54 Gd 
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v'' 






□\ i / V, 

i v □ / \\ - 







0.0 



- 160 Gd . 




D 




D /\i 







I \ 
I \ - 



-90 -60 



-30 30 
9y [degrees] 



60 90 



Fig. 4. Unperturbed angular correlations for the 2f — > Of transitions in 154 Gd 
and 156 Gd. Circles (squares) represent data taken with the detector in the posi- 
tive (negative) quadrant. The solid line is the theoretical angular correlation given 
by Eq. (4). The dotted line shows the effect of an isotropic electric- field gradient 
(Vzz = 3.44 x 10 17 V/cm 2 ). The good agreement between the data and the solid lines 
implies that the electric-field gradient is predominantly directed along the beam di- 
rection, and hence that the gadolinium foil is a quasi single crystal with the c axis 
perpendicular to the plane of the foil. 

6 Perturbed angular correlations 



6. 1 Formalism 



Since the perturbation of the angular correlation stems from changes in the 
distribution of the nuclear spins as specified by changes in the statistical ten- 
sors, the expression for the perturbed angular correlation has the same form 
as Eq. (4), with the statistical tensors (pkq) replaced by values that correspond 
to the perturbed spin distribution. The effect of the transient-field precession, 
which acts as a pure rotation around the direction of the magnetic field, is ap- 
plied first. As shown in Fig. 3, the magnetic field is directed along the y-axis. 
The statistical tensor therefore becomes 



Pkq = Y,(PkQ)D k gQ (0,A9 t{ ,0), 
Q 



(9) 
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Fig. 5. Unperturbed angular correlations for the if — > 2J transitions in 
156 ' 158 ' 160 Gd. Circles and squares represent data taken at room temperature in run 
II. The triangles represent data taken at a nominal target temperature of 90 K in 
run III. As in Fig. 4 different symbols are used to designate data taken with the 
detectors in the positive and negative quadrants. 



where (pto) are the unperturbed tensors. The sign of A# t f is reversed when 
the direction of the polarizing field is reversed. 

After the transient-field precession has been applied, a combined electric- 
quadrupole and magnetic-dipole interaction, due to the static fields in the 
gadolinium host matrix, is allowed to perturb the statistical tensor. The per- 
turbed tensors, pk> q >, are derived from the tensors, pk q , using 



(10) 



where is the perturbation factor for the combined magnetic and elec- 

tric hyperfme interactions. [G q k q k ,]* can be related to the perturbation factors 
denoted lll 9 k q k ,(P, 7), given by Alder et al. [33,34] for combined electric and 
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magnetic interactions: 



\r1l' ^(2fc+ l)(2fc'+ 1) 7T 7T QQ , k , 71 71 

i G fcfc'WJ = (27TT) ^ ^50(2' 2' ^ **' ^' 7 ^9'Q^2' 2' ^ 

(ii) 

The Euler (or spherical polar) angles 7) specify the orientation of the 
electric field gradient with respect to the magnetic field direction, which is 
along the y axis in Fig. 3. (3 is the polar angle between the directions of 
the magnetic field and the electric- field gradient. The azimuthal angle 7 is 
measured from the beam axis in the horizontal plane, i.e. the xz-plane in 
Fig. 3. Although it is not usually written explicitly, IIljj? fe ? 7) is a function 
of out and ouqt, where out is the magnetic dipole precession angle as defined 
above and ouqt is the electric quadrupole precession, where ouq is 

eQV gz , * 

W ° " 4/(2/ -l)*" (12) 



Data analysis procedures 



As a first approximation, the experimental total precession angles for the 4+ 
states, A0, can be extracted from the field-up/field-down data by conventional 
procedures [5,13] in which the experimental precession angle is related to the 
field up/down counting asymmetry e by the expression 

AG = e/S, (13) 

where S = (l/W)(dW/dO) is the logarithmic derivative of the angular corre- 
lation at the detection angle +6> 7 and 

e=(l-p)/(l+p). (14) 

The 'double ratio' p is derived from the counting rates in the detectors at ±# 7 , 
N(±0 7 ), for field up (T) and down (J.) by 



jv(+g 7 , t) JVHU) ( . 

\N{+o 1 d)N{-e 1 AY 1 0> 

Note that the factors due to integrated beam current, cross sections and de- 
tector efficiencies cancel out so that 

p = w(+e^)/w(+e^i) (is) 

and hence e is formally equivalent to 
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For the longer-lived 2j states, where the perturbations are much larger, this 
procedure does not apply. Indeed, it also underestimates the precessions of 
the 4+ states by up to ~ 10%. The procedure used here to analyze both the 
l\ and Ai data is to begin by forming the experimental double ratio, p, and 
asymmetry parameter e as usual, so that the beam current, cross section and 
efficiency factors cancel. The experimental value of wr is then extracted by 
fitting the experimental values of e (Eq. (14) and Eq. (15)) to theoretical values 
of e evaluated using Eq. (17). 



6.3 Electric quadrupole and magnetic dipole perturbations 

The combined effect of electric quadrupole and magnetic dipole interactions 
on the angular correlation and the asymmetry e is explored in Figs. 6-8. 

In Fig. 6, the transient-field precession, which is relatively small, was set to 
zero, the magnetic dipole precession angle was lot = 2 rad, and the electric 
quadrupole precession was loqt = 0.2 rad. These values are near those for 
160 Gd in the data presented below. Two extreme cases are shown. In the upper 
panels of Fig. 6 the EFG is assumed to be distributed isotropically, while in the 
lower panels it is directed along the beam axis. As discussed above, the latter 
case is very close to the real situation. Although the electric field gradient 
is directed along the beam axis, perpendicular to the initial spin orientation, 
such that it cannot perturb the angular correlation on its own, it can have 
an observable effect on the perturbed angular correlation when the magnetic 
dipole interaction is also present, because the magnetic interaction moves the 
nuclear spin out of the plane perpendicular to the beam. Thus the perturbed 
angular correlations shown in the left panels of Fig. 6 differ, depending on the 
direction and magnitude of the electric field gradient. 

Figures 7 and 8 show the dependence of the asymmetry e(65°) on the magni- 
tude of the magnetic perturbation lot for a given strength of the quadrupole 
interaction. In these figures a realistic value of the transient-field precession 
was included, A9 = —40 mrad, where the negative sign applies for 'field up'. 
The value loqt = 0.2 rad was chosen because it is slightly larger than the 
value applicable for 160 Gd, the largest considered here; loqt = 0.08 rad is 
slightly larger than the value for 154 Gd. The effect of the electric-field gradient 
is apparent in the region up to lot 1 rad when the EFG is isotropically 
distributed. However when lot » loqt, the asymmetries, e(65°), are hardly 
affected, especially when the EFG is directed along the beam direction. 

In the following analysis lot values are extracted assuming the EFG is directed 
along the beam direction, as was shown to apply for our foils in Sections 4 
and 5. loq was evaluated assuming the electric-field gradient from Ref. [2] 
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Fig. 6. Perturbed angular correlations for a hypothetical 2~[ — ► 0J~ transition which 
experiences a combined magnetic-dipole and electric-quadrupole interaction. The 
left-hand panels show the angular correlation, W(9-y), while the right-hand panels 
show the asymmetry, e(6L), derived from double ratios as described in the text 
concerning Eq.(14) and Eq.(15). In all panels the magnetic precession is lot = +2 rad 
for 'field up'. The field directions are indicated by the arrows. Solid lines show 
the case where an electric-field gradient, which causes a quadrupole precession of 
ujqt = 0.2 rad, is also present. For reference the dotted lines show the ujqt = 
case. In the upper two panels the electric-field gradient is distributed isotropically, 
while in the lower panels it is directed along the beam axis, as is the case in the 
present measurements. 

and the experimental quadrupole moments in Table 4. The parameters of the 
EFG are not critical, however. In all cases considered here essentially the same 
lot values would be extracted if EFG effects were assumed isotropic or even 
ignored altogether. 

6.4 Nuclei on damaged sites 

The analysis procedures described above implicitly assume that all of the 
implanted nuclei experience the same hyperfine magnetic field and the same 
electric- field gradient. It is well known, however, that after implantation into 
metals typically 5 - 10% of the implanted nuclei reside on damaged sites and 
therefore do not experience the same hyperfine fields as those on substitutional 
sites. 
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12 3 

cot [radians] 

Fig. 7. Asymmetry, e(65°), for a typical 2 + — ► + transition, as a function of the mag- 
netic precession angle lot when an electric-quadrupole precession of loqt = 0.2 rad 
is present, either distributed isotropically or directed along the beam direction (as 
is the case in the present work). The case of loqt = is shown for reference. In the 
present work the wr values lie between ~ 0.8 and ~ 2.2 rad, with loqt = 0.2 being 
applicable when lot ~ 2.5 rad. See also Fig. 8. 




12 3 

cot [radians] 

Fig. 8. As for Fig. 7, but for lvqt = 0.08 rad. In the present work ujqt = 0.08 rad 
is associated with lot ~ 0.8 rad. 
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Fig. 9. Asymmetry, e(65°), for typical 2 + — ► + and 4 + — ► 2 + transitions, as a 
function of the magnetic precession angle lot when a fraction of the implanted 
nuclei reside on field- free sites. Calculations are shown for field- free fractions from 
0% to 15% in 5% steps. 

It will become apparent in the following that an analysis of the present data 
assuming a unique implantation site leads to a contradiction between the 
effective fields experienced by the 4 + and 2 + states. A two-site analysis is 
necessary to resolve this contradiction. 

For the analysis of integral perturbed angular correlations it is usually suffi- 
cient to assume a two-site model in which the damaged sites have no static 
hyperfine magnetic field. This assumption will be adopted here with the ad- 
ditional assumption that there is also no net electric field gradient on the 
damaged sites. It will become apparent in the following discussion that these 
assumptions are not critical because the number of damaged sites is small. 

In a two-site model with a fraction of nuclei / on field-free sites after im- 
plantation, the perturbed angular correlation, W(A9 t f , ujt, ujqt) is replaced 



Note that the nuclei that end up on field-free sites still experience the transient 
field as they slow in the ferromagnetic medium, but they do not experience 
any hyperfine interactions after they come to rest. 

The effect of a fraction of damaged sites on the observed asymmetry e(65°) is 
demonstrated in Fig. 9. In the present fraction of nuclei on field-free 

sites, which is not taken into account, gives an apparently enhanced B st for 
the 2 + states and an apparently reduced B st for the 4 + states. The difference 
comes about because the longer-lived 2 + states are in the region where an 
increase in lot results in a decrease in the magnitude of e whereas for the 4 + 
states an increase in lot results in an increase in the magnitude of e. 



by 



(1 - /) W(Ae tf , oot, lo q t) + fW(A0 t{ , 0, 0). 
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Table 7: Summary of 2 + precessions. 



Run 


A0 tf 


{jJ Q T 


e(65°) x 10 3 


ut(6 


5°) e 


(120°) x 10 3 




120°) 


H 




(mrad) 


(mrad) 




(mrad) 




(mrad) 


(mrad) 










1-site 


2-site 




1-site 


2-site 


1-site 


2-site 












154 Gd 












II 


-35 


67 


-204(14) 


904l?| 


69811 


221(16) 


7841H 


7+90 

°° ' -103 


8491H 


640±1| 


III 


-43 


67 


-193(10) 


95611 


828l^ 7 


180(16) 


10301^1 


865llJ 3 


972lg 


8361^ 


IV 


-43 


67 


-187(11) 


995l 8 7 ° 


8631^ 


153(13) 


124211H 


10731H 


10621*;° 


9181S 












156 Gd 












II 


-31 


139 


-148(8) 


12991% 


10391^ 


130(7) 


1529!^ 


11681^ 


139711 


109811? 


III 


-39 


139 


-120(5) 


1635l 8 7 ° 


14351! 


117(7) 


170611^ 


1456l^ 5 


165711 


144ll| 


IV 


-39 


139 


-107(5) 


18581^° 


16261^ 


113(7) 


177511 33 


151411^ 


18271^ 


158411 












158 Gd 












II 


-31 


165 


-123(11) 


15881H 


1274liS 


107(8) 


18931H 


14471H 


17421H 


136711 


III 


-38 


165 


-100(6) 


19841H 


1740±1S 


94(9) 


21681^ 


1848l? 7 ° 


20291H 


1767l^ 4 


IV 


-38 


165 


-91(8) 


21981?! 


19221^ 


101(8) 


20011H 


170911^ 


208411H 


17971H 



Table 7: (Continued) 



Run A0 tf u q t e(65°) x 10 3 cur(65°) e(120°) x 10 3 cjt(120°) (lot) 

(mrad) (mrad) (mrad) (mrad) (mrad) 

1-site 2-site 1-site 2-site 1-site 2-site 



160 Gd 



II 


-30 


177 


-111(12) 


1770l?| 


1424ll| 


101(10) 


20101^ 


15411H 


1885llJo 


1483±1^ 


III 


-36 


177 


-87(7) 


2298±?Si 


20161H 


112(11) 


1774l?| 


15241H 


203611^ 


1768ll?J 


IV 


-36 


177 


-73(8) 


277011° 


2415111! 


104(10) 


19291^1 


16561H 


217911^ 


18851H 




Fig. 10. Perturbed angular correlations for the 2^ — > 0^ transitions in Gd and 
156 Gd measured in run II. The solid curves are calculated for the average single-site 
lot values in Table 7. These data make it clear that the ojt values for the 2f states 
are much larger than 0.5 rad. 

6.5 Results: 2t states 



Table 7 summarizes the results of the present measurements of the static- 
field precessions for the 2f states in i 54 > 1 56,158,i60q ( j_ rpj^ ^ a ^ a were analyzed 
adopting both a single-site model (i.e. / = in Eq. (18)) and a two site model 
(/ 7^ 0). The procedure by which the field-free fraction was determined for 
the two site-model analysis is described in section 6.7. 

While the ujt values Table 7 were extracted exclusively from the asymmetries, 
e(65°) and e(120°), these values are consistent with the full perturbed angular 
correlations, where measured. For example, the perturbed angular correlation 
data for 154 Gd and 156 Gd measured in run II are shown in Fig. 10, and com- 
pared with the perturbed angular correlations corresponding to the ujt values 
for the single-site model (see Table 7). 

To some extent the extracted magnetic-dipole precession angles, ujt, depend 
on the assumed transient-field precession angle, A9 t{ , and the quadrupole pre- 
cession, ujqt. Figures 6 and 8 discussed in the previous subsection show that 
ujt is very insensitive to reasonable assumed values of ujqt. The effect of the 
transient-field contribution is somewhat counter intuitive at first sight: ujt 
decreases by ~ 20 mrad when the magnitude of A0 t f (which is negative) in- 
creases by 10 mrad. (See Fig. 4 in Ref. [35] for a plot of e versus ojt showing 
the effect of the transient-field contribution.) Because the ujt values are so 
large, however, an increase in A# t f by ~ 50% would lead to a change in ujt 
for 154 Gd ( 160 Gd) by only ~ 4% (~ 2%). 

A fixed value of A9 t f/g = —100 mrad was adopted for run IV, and scaled 
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according to the foil magnetization for runs II and III as obtained in a prelim- 
inary fit to the data. This value was chosen because it is consistent with the 
previous data [1] once a correction is made for the difference in Gd recoil ve- 
locities due to the difference in the beam energies. In principle, A9 t { could be 
obtained from the data for the 4 + states. Unfortunately, however, the present 
data for the 4 + states are not sufficiently precise to fit both the static- and 
transient-field precessions as free parameters (see below). 

Results presented for the two-site model assume a field- free fraction of 11.6% 
in run II and 6.4% in runs III and IV. These values were determined by 
requiring consistency between the extracted static-field strengths for the 2f 
and 4f states, as described in section 6.7. 

The effective static-field strengths derived from the average precession angles 
given in Table 7 are summarized in Table 8. As will be discussed below, the dif- 
ferences in these effective static-field strengths from run to run are attributed 
to changes in the magnetization of the gadolinium foil target associated with 
different foil textures and different beam heating effects. 
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Table 8 

Effective static-field strengths derived from 2 + -state precessions. 



Nucleus B st (Tesla) Single-site analysis B st (Tesla) Two-site analysis 

Run II (A) Run III (A) Run IV (B) Run II (A) Run III (A) Run IV (B) 

154 Gd -24.1(2.4) -27.6(2.5) -30.2(2.8) -18.2(2.1) -23.8(2.2) -26.1(2.5) 

156 Gd -23.1(1.1) -27.4(1.2) -30.2(1.4) -18.1(0.9) -23.8(1.0) -26.2(1.2) 

158 Gd -25.6(1.8) -29.8(1.8) -30.6(2.1) -20.1(1.4) -26.0(1.5) -26.4(1.8) 

160 Gd -27.7(2.7) -29.9(2.6) -32.0(3.1) -21.8(2.1) -26.0(2.2) -27.7(2.6) 

average -24.3(0.8) -28.3(0.9) -30.5(1.0) -19.0(0.7) -24.6(0.8) -26.4(0.9) 
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Fig. 11. Left: Previous precession results, AQ/g, for the and 6^ states in the 
even Gd isotopes [1] plotted as a function of the mean level lifetime r. These data 
were obtained following Coulomb excitation with 56 MeV 16 beams. Right: Present 
data for 4+ states in i56,i58,i60 Gd from 

run IV, obtained with 40 MeV 16 beams. 
In both panels, the solid line, for B st = —25 T, is the best fit to the data. The 
dotted line forces the static field strength to be —30 T as implied by the single-site 
analysis of the data for the 2f states. 



6.6 Results: 4t states 



Table 9 summarizes the results of the present measurements on the states 
in 156 > 158 > 160 Gd. The 4 + state in 154 Gd was populated too weakly in the present 
work to be analyzed. The precession angles ujt were obtained by the rigorous 
procedure described in section 6.2 assuming the same A0 t f values as adopted 
for the analysis of the 2± states (Table 7). 

If both sides of Eq. (1) are divided by g it becomes 

AQ/g = (u/g)T + A6 t{ /g. (19) 

From this equation it can be seen that a plot of AQ/g values versus r should 
lie on a straight line with a slope u/g, which is proportional to B st , and with 
an intercept A6 t f/g at r = 0. To aid comparison with previous work, the final 
column of Table 9 shows the total precessions, A©, according to the single- 
site model. The left panel of Fig. 11 shows the previous data, from [1], for 
the 4+ and 6^ states in the even Gd isotopes so plotted. For comparison the 
right panel shows the present data for the 4+ states from run IV (single-site 
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Table 9 



Summary of 4f 


state 


precessions. 










Run 


Target 


A% 


e(65°) 


e(120°) 


H 


{AG)/g 








(xlO 3 ) 


(xlO 3 ) 


(mrad) 


(mrad) 












1-site 


2-site 


1-site 










156 Gd 








II 


A 


-31 


-40(12) 


1(12) 


50(9) 


56(10) 


52(25) 


III 


A 


-39 


-21(10) 


24(14) 


61(8) 


65(9) 


58(21) 


IV 


B 


-39 


-26(10) 


5(12) 


56(8) 


60(8) 


44(21) 










158 Gd 








II 


A 


-31 


-38(10) 


35(10) 


68(7) 


77(9) 


102(20) 


III 


A 


-38 


-55(8) 


59(12) 


97(8) 


105(9) 


155(21) 


IV 


B 


-38 


-65(8) 


65(10) 


108(8) 


117(9) 


181(21) 










160 Gd a 








II 


A 


-30 


-36(10) 


49(11) 


72(8) 


82(9) 


115(23) 


III 


A 


-36 


-64(9) 


48(13) 


96(9) 


104(10) 


157(24) 


IV 


B 


-36 


-84(9) 


64(10) 


117(9) 


126(10) 


209(23) 



a A correction has been applied to account for a ~ 8% contribution to the 
intensity of the 4+ ->■ 2+ transition in 160 Gd due to the 9/2~ ->■ 7/2~ 
transition in 157 Gd. 

analysis). The results of the previous work [1] have been reanalyzed using the 
present adopted g factor and lifetime values (Tables 3 and 4). These previous 
data were fitted treating both ui jg and A9 t f/g as free parameters. The result, 
shown as the solid line in the left panel of Fig. 11, corresponds to B st = 
—25(3) T, in agreement with the value of —25 T obtained in run IV (single- 
site analysis). 

Unfortunately the present data for the 4f states are not sufficiently precise 
to fit both the slope and the intercept as free parameters. The transient-field 
precession was therefore set to the values adopted in the analysis of the 2 + 
states (see Table 7). These values and the extracted static- field strengths for 
the one- and two-site models are presented in Table 10. 

Since the previous measurements [1] were made with 56 MeV 16 beams, 
somewhat higher than the 40 MeV beams used here, the transient-field con- 
tribution must be larger in the previous work. According to the Rutgers [36] 
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Table 10 

Hyperfine fields extracted from 4+ state precessions. 



Run Target A0 ti /g B st 
(mrad) (Tesla) 









1-site 


2-site 


Ref. [1] 




-136(26) 


-25(3) 




II 


A 


-80 


-16.7(1.3) a 


-19.0(1.4) 


III 


A 


-93 


-22.7(1.4) a 


-24.4(1.5) 


IV 


B 


-100 


-24.7(1.4) a 


-26.7(1.5) 



a The uncertainty in the transient-field contribution is not 
included in the quoted error. 

and Chalk River [6] parametrizations, the expected difference in A6 t {/g is 
about 30 mrad. As shown in Table 10, the present data for target B in run 
IV (single-site analysis) are consistent both with the expected difference in 
transient-field precession and with the same static field strength as observed 
in Ref. [1] . The present measurements from run IV are also in agreement with 
the static- field strengths obtained by Hausser et al. [6]. 

All of the previous work has assumed a single-site for the implanted nuclei. 
However a comparison of the results in Tables 8 and 10 shows that the as- 
sumption of a single implantation site implies that the effective static-field 
strengths for the 4+ states are significantly smaller than those experienced by 
the 2t states. 



6. 7 The field-free fraction 

Figure 12 illustrates the procedure used to determine the field-free fraction /, 
according to the two-site model described in section 6.4. In order to achieve 
consistent B st values for the 2 + and 4 + states requires that in run II, 11.6^3% 
of the implanted nuclei reside on field-free sites. In runs III and IV, the field- 
free fractions are 6.Tt\'^,% and 6.lt\'%%, respectively, giving an average of 
6.4 ± 1.3%. Thus the two-site analysis adopted / = 0.116 for run III and 
/ = 0.064 for runs III and IV. The initial expectation would be that / should 
be the same for all three runs. It is not clear why a higher percentage of 
ions apparently reside on damaged sites in run II. Since the hyperfine fields 
experienced by the implanted nuclei in run II are very different from those 
in runs III and IV, this difference might indicate that there is another effect 
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Fig. 12. Determination of the field-free fraction for run IV. 



associated with the perturbed angular correlations that has not been taken 
into account. The apparently different value of / was retained for the two-site 
analysis of run II, subject to the caveat that it might not represent the true 
field-free fraction. 



7 Discussion 

7. 1 Synopsis 

The primary motivation for the present measurements was to obtain an in- 
beam measure of the local magnetization of the gadolinium target foils, which 
can vary from foil to foil and for different beam-heating conditions. Clearly, 
the effective hyperfine field strength varies with the magnetization. The fol- 
lowing discussion explores the extent to which ratios of effective fields can be 
interpreted as ratios of the host magnetizations. 

Along with the presentation of the experimental data, the previous sections 
have established that the gadolinium foils used in the present experiments 
are quasi single crystals in which the electric-field gradient is directed along 
the beam direction (sections 4 and 5.2). It has been noted that the effect of 
the electric-field gradient is negligible for 4f states. It has also been shown 
that, for the analysis procedures adopted here, the effect of the electric-field 



29 




cot [radians] cox [radians] 



Fig. 13. Asymmetry, e(65°), as a function of lot showing the effect of misalignment 
between the external field direction and the internal hyperfine field direction. The 
angles (3 and 7, which specify the direction of the internal field, are defined as in 
Fig. 3. For the 4 + —* 2 + case the curves for = 30° on the cone and 7 = 0° are 
almost indistinguishable. 

gradient along the beam direction is essentially negligible for the longer-lived 
2 1 states as well. It follows that the presence of the electric-field gradient 
does not impede the use of the effective static-field strength as a measure of 
the local magnetization. 

Before coming to a discussion of the results in terms of the local magnetization 
and beam heating effects (section 7.3), it is necessary first to discuss the effect 
of the hyperfine field being misaligned with respect to the external field due 
to domain misalignment below the saturation magnetization (section 7.2). 

7.2 Domain rotation below the saturation magnetization 

It has been demonstrated for several impurities implanted into iron that inter- 
nal hyperfine magnetic fields may be misaligned with respect to the direction 
of the external polarizing field, and that this effect is associated with domain 
rotation in an incompletely saturated sample [37]. At some level the effect is 
expected for all impurity-host combinations. 

Figure 13 shows the effect of misaligned hyperfine fields on the analysis of 
the perturbed angular correlations. Two limiting cases of misaligned fields are 
assumed. In the first case the internal fields are assumed to lie on a cone of 
half angle f3 = 30°, whose axis is the external field direction. In the second 
case the internal field is assumed to lie in the plane denned by the beam axis 
and the external field direction, i.e. at angles j3 = 30°, 7 = 0° in terms of the 
co-ordinate frame in Fig. 3. Given that the gadolinium target has a texture 
such that it resembles a quasi single crystal with the c axis along the beam, 
the latter case is likely to be more realistic in the present work. The formulae 
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for evaluating these perturbed angular correlations are given in the Appendix. 

It can be seen from Fig. 13 that if the internal field is misaligned with respect 
to the external field, then the true precession angle around the internal-field 
direction is always larger than that derived when no misalignment is assumed. 

At least part of the difference between the effective static fields observed here 
and the Mossbauer result at 4 K (B st = 37.3(5) T [2]) is likely to be associ- 
ated with misalignment between the internal and external fields. Indeed the 
subtle overall tendency for the effective static fields in Table 8 to increase 
slightly from 154 Gd to 160 Gd, as the 2f state lifetimes and hence the lot val- 
ues increase, might be associated with misalignment between the internal and 
external fields, which has not been included in the analysis. (Note that the 
difference between the different curves in Fig. 13 decreases as uot increases.) 
However the trend is below the statistical precision of the data. 

It can be concluded that although the perturbed angular correlations may 
show some sensitivity to the direction of the internal field, the ratios of out 
values derived assuming that the internal and external fields are parallel can 
still be interpreted as magnetization ratios, to a very good approximation. 

7.3 Beam heating and relative magnetizations of foils 

There is no significant difference between the ratios of effective fields deter- 
mined from the one-site or two-site analysis of the data. The following discus- 
sion will therefore use the ratios of the observed precession angles in the sim- 
pler single-site analysis as the measure of the relative magnetizations. Table 11 
summarizes the results of these in-beam relative magnetization measurements. 

The only difference between runs II and III is the beam intensity, which is 
nearly an order of magnitude larger in run II (see Table 1). If it is assumed 
that the local temperature at the implanted nuclei is near 100 K in run III, 
the reduction in magnetization in run II, to ~ 85% of the value near 100 K, 
must correspond to a significantly higher local temperature of ~ 190 K (see 
Fig. 2). 

Thermal conductivity calculations were performed with the QuickField pack- 
age [38] to investigate the effects of beam heating on gadolinium target foils, 
with and without copper backing layers. Bulk thermal conductivities were as- 
sumed for the gadolinium and copper layers. Some results are presented in 
Fig. 14 and Fig. 15. The temperature in the beam spot according to these 
calculations is included in Table 11. 

Despite the schematic nature of these calculations, and the fact that the foils 
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Table 11 

Relative magnetizations from in-beam precession measurements. 



Kun, larget 




D B a 




Relative 


magnetization 




(W) 


(10 9 W/m 3 ) 


(K) 


4+ 


2 + average 


II, A (Gd) 


0.1 


1.5 


~ 185 


0.57(8) 


0.80(3) 0.77(3) 


III, A (Gd) 


0.015 


0.2 


~ 90 


0.83(9) 


0.93(4) 0.91(4) 


IV, B (Gd + Cu) 


0.1 


2.3 


~ 90 


1 


1 1 



a Pb is the power deposited in the target by the beam. Db is the power 
density in the beam spot. 

b Calculated temperature at the center of the beam spot. 

c Relative magnetizations determined from ratios of precession angles given 

in Tables 7 and 9. 




20 urn Gd 1 jim Gd + 5 urn Cu 



Fig. 14. Schematic thermal conductivity calculations. The front view of the target 
is shown on the left. To the right are shown side views of the targets which indicate 
the temperature profiles between the center of the target and the point of contact 
with the target frame. The two cases shown approximate target A in run II (left) 
and target B in run IV (right). 

are rather thin (so bulk thermal conductivities may not apply), the calcula- 
tions are in very good qualitative agreement with experiment. The two exam- 
ples shown in Fig. 14 were chosen to resemble target A in run II and target B 
in run IV. The schematic calculations reproduce the difference in temperature 
required by the experimental data. Fig 15 shows the linear variation in temper- 
ature for a 20 /zm thick gadolinium foil as the beam power is changed. Again, 
these schematic calculations reproduce the experimental difference between 
runs II and III. 

There remains a difference, of the order of 10%, between the magnetizations of 
target A and target B under conditions where the beam heating is essentially 
negligible (runs III and IV). This difference must be attributed to differences 
in the crystalline texture of the foils. 

The thermal conductivity calculations also suggest that even a relatively thin 
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200 




50 1 1 1 1 1 

1 2 

Power density (10 9 W/m 3 ) 

Fig. 15. Calculated temperature at the beam spot as a function of beam 
power-density, Db, for a 20 fim thick gadolinium foil. 

layer of copper evaporated onto a gadolinium foil will greatly assist in keeping 
the temperature in the beam spot near liquid nitrogen temperature, under 
conditions where the beam intensity cannot be kept very low. Such a procedure 
was used for this purpose in Ref. [10]. 



8 Summary and Conclusions 

The effective hyperfme fields experienced by Gd ions recoil-implanted into 
gadolinium foils have been measured for two targets (one copper-backed) and 
with differing beam intensities. The effects of (i) the transient-field interac- 
tion, (ii) electric-field gradients, and (iii) nuclei residing on damaged, field-free 
sites have been evaluated. The possible effects of a misalignment between the 
external field and the direction of the internal magnetization was discussed. It 
was found that the effective hyperfme magnetic field varies from target to tar- 
get and with the power deposited by the beam, particularly when the target 
is an un-backed gadolinium foil. 

To a very good approximation, the ratios of hyperfme magnetic fields can be 
interpreted as ratios of the host magnetization. The changes in hyperfme field 
strength, and hence magnetization, with beam intensity can be correlated with 
the expected temperature rise in the beam spot due to the power deposited by 
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the beam. A layer of copper evaporated onto the gadolinium foil can greatly 
enhance the dissipation of beam power and minimize this temperature rise. 

The results of the present measurements were used in a recent study of the 
transient-fields for high- velocity Ne and Mg ions traversing gadolinium [16], to 
correct for differences in the magnetizations of the gadolinium foils in targets 
A and B. 
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Appendix 



In this appendix the formulae are given for calculating perturbed angular 
correlations in the case where the internal hyperfine fields are misaligned with 
the external field. Since a full quantitative analysis is not required here, it is 
necessary to consider the case where only q = terms are non-zero in Eq. (4). 
It is then possible to define 

a k = V2k + l( Pk0 )F k Q k . (20) 



The case where the internal field is assumed to be distributed equally on a 
cone of half-angle (3 to the external field has been considered previously by 
Ben-Zvi et al. [39]. Their expression can be rewritten in the form 



W{6 1 ,P) = Y j a k [d^)} 2 \d% q m 2 

kQq 

cos(gfl 7 - Q(A6 Q + Afl tf )) 

where 

tan(QA6> Q ) = Qcut, (22) 

and for convenience in the present application, the spherical harmonics in 
their expression have been replaced by the dq q {P) matrices for the second 
Euler rotation (about the y axis). 

To our knowledge, the more general case, where the internal field has a single 
specific orientation to the external field, specified by the spherical polar co- 
ordinates 7) (see Fig. 3), has not been given previously. The result is 



fc<2<?l<32 

cos(g 2 fl 7 - - g l ) 1 - Q(A6 Q + A0 tf )) 

This expression reduces to Eq. (21) when 7 is averaged between and 2ir 
radians: 

1 r 2lT 

W(6 y ,/3) = — W(0 7 ,/3,7)d7 (24) 
Ztt jo 
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